Breakdown of the power-law decay prediction of the heat current correlation in 
one-dimensional momentum conserving lattices 
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We show that the asymmetric inter-particle interactions can induce rapid decay of the heat current 
correlation in one-dimensional momentum conserving lattices. When the asymmetry degree is ap- 
propriate, even exponential decay may arise. This fact suggests that the power-law decay predicted 
by the hydrodynamics may not be applied to the lattices with asymmetric inter-particle interactions, 
and as a result, the Green- Kubo formula may instead lead to a convergent heat conductivity in the 
thermodynamic limit. The mechanism of the rapid decay is traced back to the fact that the heat 
current has to drive a mass current additionally in the presence of the asymmetric inter-particle 
interactions. 
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Onsager's classic work \l\ has shown that in the linear 
response region, the current Ji of a physical quantity can 
be expanded in terms of the so-called thermodynamic 
forces, i.e., Ji = i^ijPj- The Onsager coefficient 
which describes the coupling between the forces Fi and 
Fj , can be calculated in terms of the correlation function 
of the spontatous current fluctuations as [2j , 



iij — lim lim — 



(1) 



where C^J{t) = {(Jtit) - (J,))(J,(0) - (J,-))} and L is the 
linear dimension of the system along which the current 
flows. The celebrated Onsager reciprocal relation formu- 
lates the relation between the coupling coefficients, which 
gives Kij — Kji. However, we recall that under what con- 
dition the coupling coefficients are non-vanishing, and 
what a role the couplings may play in characterizing the 
transport process, have not been clarified. 

Equation (1) is the well-known Green-Kubo formula. 
It provides a way for calculating the transport coefficient 
by considering the current fluctuations in equilibrium 
state. The traditional hydrodynamic approach assumed 
that the current correlation decays rapidly, i.e., exponen- 
tially fast, so to ensure the convergence of the integral in 
the Green-kubo formula 0, However, after Alder nu- 
merically evidenced the 'long time tail' of the correlation 
function of the energy current in 1970 [HI in a gas model, 
a lot of theoretical analysis as well as numerical simu- 
lations have shown that the autocorrelation functions of 
currents in one-dimensional (ID) fluids may generally de- 
cay in a power-law manner instead [6- 8] . The power-law 
decay is explained within the framework of hydrodynam- 
ics, where the slow diffusion of long wave hydrodynamic 
modes or the ring-collision of particles are ascribed to be 
the underlying mechanisms. 

In recent decades, low-dimensional materials such as 
nanowires and graphene flakes have come into focus in 
many disciplines. The studies of them are undergoing 
rapid progress for both fundamental research interests 



and various intriguing applications. The heat transport 
properties of low-dimensional lattice models have been a 
particularly interesting issue. An important progress is 
that the hydrodynamic analysis has been extended to the 
study of the transport behavior in lattices. Based on in- 
tensive theoretical analysis [1, 0- 14 1 and numerical stud- 
ies (see for example 0,1a] and references there in), for ID 
momentum conserving fluids and lattices, at present it is 
generally believed that the current correlation should de- 
cay in the power-law manner. An important consequence 
of the power-law decay is that the integral in the Green- 
Kubo formula will diverge. For a flnite system size, to 
truncate the integration at a reasonable time may lead 
to a flnite transport coefficient, but, however, it may di- 
verge in the thermodynamical limit in a way of Hij ~ L'^ 
(with 1/ > 0). 

Meanwhile, several counterexamples to the divergent 
heat conductivity have also been reported, including the 
rotator model [16|, a ID lattice in an effective magnetic 
field [l3|, the variant ding-a-ling model 18 1, and lattice 
models with asymmetric inter-particle interactions [l9| . 
We notice that in the counterexamples studied in Ref. 
[l8| and [l^, asymmetric inter-particle interactions are 
quite relevant. Hence a natural question is what effect 
the asymmetric interactions may have on the 'long time 
tail' decay and result in a finite, convergent heat conduc- 
tivity. In fact, lattice models with asymmetric interac- 
tions have been studied in the literature, both analyti- 
cally 0, M m and numerically [lol. [23] . In particular, 
the hydrodynamic analysis based on the Burgers equa- 
tion suggests that a ~ 1/2 for systems with asymmetric 
interactions and a = 1/3 for those with only symmet- 
ric interactions Q , which agrees with the result based on 
the Zwangzig-Mori equation and the self-consistent mode 
coupling theory [l3] ■ 

In this paper, by employing two paradigmatic lattice 
models, the Fermi-Pasta-Ulam-Q-/3 (FPU-a-/3) model 
and the Lennard-Jones (L-J) model, we show that the 
power-law decay of the current correlation may vanish 
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FIG. 1: Plots of the inter-particle potentials investigated: The 
potential of (a) the FPU-a-/3 lattice model and (b) the L-J 
lattice model. 



in the presence of the asymmetric interactions. Instead, 
rapid decay fast than the power law is observed in both 
models. It is well known that systems with symmetric 
and asymmetric interactions have remarkable different 
properties. For example, the thermal expansion effect is 
exclusive in the latter. In the following, we shall conjec- 
ture that the thermal expansion effect is also responsi- 
ble for the rapid decay of current correlation functions. 
For this aim we have to confront the previous numerical 
results in the literature and the contradictory hydrody- 
namic predictions. On one hand, we shall provide careful, 
large-scale simulation results, and show the importance 
to keep the lattice features of the models in the simula- 
tions and the sensitive dependence of the results on the 
asymmetry degree. As to the hydrodynamic predictions, 
we think it is still open. As the hydrodynamic analy- 
sis is based on the linearized hydrodynamic equations, 
one possibility is that the role of the asymmetric inter- 
actions may not have been sufhciently considered if their 
main effects are contained in the high order terms. We 
shall show that the energy current can excite the mass 
current in systems with asymmetric interactions, while 
there lacks such a coupling between the energy and the 
mass current in systems of symmetric interactions. We 
shall also illustrate that, different from the fluid models, 
the lattice structure can additionally scatter the currents. 
These facts are crucial for the rapid decay of current cor- 
relation but have not been involved explicitly in the pre- 
vious hydrodynamic analysis. 

Our models are defined by the Hamiltonian 
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where pi and Xi are the momentum and position of the zth 
particle, respectively, and V the potential between neigh- 
boring particles. The component particles are assumed 
to be identical and have unit mass, and the lattice con- 
stant is set to be unity so that the system length L equals 
the particle number N. The inter-particle interactions in 



1 9 a T 1/1 
V{x) - -x^ - -x' + -x\ 



(3) 



where the parameter a controls the degree of asymme- 
try as illustrated in Fig. 1(a). For a = the system 
reduces to the Fermi-Pasta-Ulam-/3 (FPU-/3) model with 
symmetric potential only. To well reveal the effects of 
the asymmetry, in our simulations the average energy 
per particle, denoted by e, is fixed to be e = 0.1 such 
that the averaged potential energy per particle is about 
0.05. 

The potential of the L-J model is 



V{x) = [(- 



1 



1 



1 



_ 2( )" + 1] 
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(4) 



This potential may well approximate the inter-particle 
interactions in many real materials, and hence has impor- 
tant practical implications. In this potential the degree 
of asymmetry is controlled by the parameter set (m, n). 
(See Fig. 1(b) for several examples.) It should be noted 
that if the average potential per particle is larger than 
one then the potential models fluids instead. So in our 
simulations the average energy per particle is fixed to be 
£ = 0.5 to make sure that our model is a lattice. As 
shown in Fig. 1(b), the degree of asymmetry increases as 
{m,n) changes from (12^6) to (6,3), and to (2, 1). 

dv 



The energy current [21"] is defined as = i. 
For a lattice the energy current is equal to the heat cur- 
rent because there is no residual global velocity 01 ■ To 
numerically measure the current in the equilibrium state, 
the system is first evolved from an appropriately assigned 
random initial condition for a long enough time (> 10*) 
to ensure that it has relaxed to the equilibrium state; 
then the current at ensuing times is measured. The pe- 
riodic boundary condition is applied in the calculations. 

Figure 2 shows the autocorrelation of the energy cur- 
rent Cjqjq{t) = Ij^IoIj^Io)) for the two models. The re- 
sults are presented in semi-log scale in panel (a) and (c) 
and in log-log scale in panel (b) and (d). In generating 
Fig. 2, the system size is fixed to be = 2048 which is 
sufficiently large because we have checked that the results 
does not change if the system size is increased further. It 
can be seen that in the FPU-a-/3 model which involves 
symmetric potential, the correlation function decays in a 
power law Cj j it) ^ with 7 ^ 0.67, which agrees 



well with previous studies [21| . Therefore, it provides an 



example that with the symmetric interactions, the long- 
time-tail prediction applies. However, for the FPU-a-/? 
model with a > 1, we can see that the decay become 
faster than the power-law manner which can be roughly 
regarded to be exponentially. For the L-J model, the 
decay is much close to a perfect exponential manner for 
different sets of control parameters (m, n). With the in- 
crease of the asymmetry, the decay become faster and 
faster. 
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FIG. 2: The autocorrelation function of the heat flux, (a) 
and (b) are respectively the semi-log plot and log-log plot of 
FPU-a-/3 model; (c) and (d) are the same but for the L-J 
model. 



The thermal expansion effect induced by the asym- 
metric interactions could be a key for understanding the 
mechanism (IS] . Thermal expansion implies the coupling 
between the energy and mass distribution. Because of the 
coupling, redistribution of energy can induce redistribu- 
tion of mass. In the following, we demonstrate that in 
equilibrium state, asymmetry interactions can do result 
in the coupling between the fluctuations of the energy 
current and the mass current. As the global mass cur- 
rent in equilibrium state is zero due to the fact that the 
total momentum of the system is zero, the cross correla- 
tion between the global energy and mass currents is zero 
as well. However, the local mass current of a part of the 
system, i.e., Jm = Mv, may still fluctuate. Here M is the 
mass of the whole part of the system and v represents its 
center-of-mass velocity. We thus investigate the coupling 
between the local energy current and the local mass cur- 
rent. To be concrete, in the simulations we consider a 
quarter of the whole system as our local subsystem. 

Figure 3(a) shows the cross correlation of the local en- 
ergy current Jq and the mass current Jm as function of 
time for the FPU-a-/? model with a — 2. It can be 
seen that {Jq{0)Jm{t)} and {Jm{t)Jq{0)) perfectly agree 
with each other, suggesting that the Onsager reciprocal 
relation, Kqm — Kmq, holds even for the subsystem. In 
Fig. 3(b), the cross correlation at the same time, i.e., 
{Jq{t)Jra{t)), IS showu as a function of the asymmetry 
parameter a. It can be seen that the cross correlation 
decreases with the decrease of the asymmetry. In partic- 
ular, in the symmetric limit with a = 0, the correlation 
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FIG. 3: (a) The cross correlation functions between the lo- 
cal heat flux Jq and the local mass flux Jm for the FPU-Q-/3 
model with a = 2. (b) The dependence of the cross correla- 
tion function {Jg(0)Jm(0)) on the asymmetry parameter a. 

vanishes, which reveals clearly that it is the asymmetry 
that sustains the coupling between the currents. 

The fact that a symmetric potential leads to a vanish- 
ing cross current correlation can be shown analytically. 
In the continuum limit, the local energy current can be 
written as j = {i^}- In addition, the mass density is a 
scalar variable and satisfies p{x) = p{—x) because the po- 
tential is symmetric and hence we have V{x) = V{—x). 
Considering all of these, it is straightforward to derive 
that 

{j x ^^J-^^ xpix)dx) = - {J x^^^xp{x)dx), 

which implies a zero cross correlation immediately. Here 
the integration is taken over the subsystem on which the 
local mass current is defined. Here (—1,1) represents the 
lattice segment of the subsystem. 

Therefore, in a system with the asymmetric interac- 
tions, once a fiuctuation of energy current forms sponta- 
neously, a mass density fluctuation will be excited. The 
resultant mass current is driven by the energy current 
due to their nonzero cross correlation, hence the mass 
current should serve as a resistant factor to the energy 
current. 

This mechanism alone, however, can not explain the 
rapid decay of the current correlation. The conversional 
hydrodynamic theory has predict that low-dimensional 
fluids should show slow-decay of current correlation. For 
the ID diatomic gas previous simulations have shown 
the power-law decay of Cqq{t) ^ t^^-^"^ . We have reper- 
formed the simulation and obtain the same result. In 
this model, particles interact via the hard-core collisions, 
which is the limit case of the asymmetry potential. This 
fact indicates that the lattice feature should also be cru- 
cial to the fast decay. 

The spatiotemporal cross correlation of fluctuations of 
the local energy current and the local mass current may 
shed light on the role of the lattice feature. For this 
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FIG. 4: The spatiotemporal cross correlation of fluctuations 
of the local energy current and mass current for (a) the FPU- 
a-P model (Green) and (b) the gas model (Blue). 



purpose we divide the whole system into N equal bins. 
Then the global current J can be expressed as the sum 
of the local ones, i.e., J = J2k jk, where jk represents the 
local current of the A:th bin. We use 



C{x,t) 



it)) 



(5) 



to measure the spatiotemporal cross correlation between 
the local energy current jl in the fcth bin at t = and 
the local mass current j[" in the Ith bin at time t, where 
X = {I — k). Figure 4 shows the results for the FPU- 
a-(3 model with a = 2 and e — 0.1 and the diatomic 
gas model with mass ratio 1/3 and and e = 1. The dif- 
ference is remarkable: In the lattice model the positive 
correlation peaks are followed by regions with negative 
correlation, while in the gas model there are only posi- 
tive correlation peaks. These observations indicate that 
in the lattice model, as the mass current flows in one 
direction, a companying mass current in the reverse di- 
rection will be excited as well. Because of the coupling 
between the energy and mass current, it also implies that 
an energy current in the reverse direction is induced as 
well. In other words, the currents are reflected as they 
flow forwards, and a decay mechanism is thus resulted in. 
However, such a mechanism is absent in the gas model. 

In summary, we have shown that the asymmetric inter- 
particle interactions could result in nonvanishing cou- 
pling between the local energy current and the local mass 
current, and the coupling coefficients follow the Onsager 
reciprocal relation. In the case of the symmetric inter- 
actions, the coupling vanishes. The asymmetric interac- 
tions can induce rapid decay of the current correlation in 
ID momentum conserving lattices, and results in finite 
thermal conductivities in the thermodynamic limit. As 
the asymmetric interactions, signalized by the thermal 



expansion effect 25|, are common in real materials, this 
finding should have universal implications. 

Our simulations indicate that even perfect exponen- 
tial decay can be observed with appropriate degree of 
asymmetry. However, based on numerical studies with 
finite-size effects, at present it is difficult to conclude with 
certainty that the asymmetric interactions can generally 
lead to the faster than the power-law decay in the cur- 
rent correlation. The degree of asymmetry is not only 
determined by the control parameters but also depends 
on the temperature of the system. For example, at the 
same control parameter of a = 1 we have observed fast 
decay of the current correlation in the FPU-a-/3 model, 
while in previous studies [3] a power-law like decay was 
observed. The difference roots in the temperature: in 24 1 
a much higher temperature was investigated, hence the 
potential is in effect dominated by the quartic term and 
shows a more symmetric structure. To clarify in detail 
the dependence of the decay behavior of current corre- 
lation on the degree of asymmetry should be a task for 
next studies. 

Another more challenging task for future studies is to 
answer why the direct simulation results as presented 
here disagree with certain hydrodynamic predictions. 
One clue we provide in this paper is that the coupling 
effect between the heat and the mass current. With the 
asymmetric interactions, local heat currents may excite 
local mass currents, hence the latter may be an resistant 
factor to the former. In the case of symmetric interac- 
tions, the coupling effect vanishes. On the other hand, 
the currents can also be scattered by the lattice struc- 
ture in the presence of asymmetric interactions. This 
mechanism is absent in the fiuid model we studied. This 
observation may explain why in the one-dimensional di- 
atomic gas a power-law decay can be observed though 
the inter-particle interactions are also asymmetric. 
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